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SUMMARY 


A comparison was made be*:ween experimental pressure distributions measured during 
testing of the Vought A-7 fighter and the theoretical predictions of four transonic 
potential flow codes. Isolated wing and three wlng-body codes were used for compari- 
son. All comparisons are for transonic Mach numbers and include both attached and 
separate flows. 

In general, the wing-body codes gave better agreement with the eifperlment than 
did the Isolated wing code but, because of the greater complexity of the geometry, 
were found to be considerably more expensive and less reliable. 


INTRODUCTION 


The past 10 years has seen the development of numerous transonic-flow analysis 
codes based on finite difference solutions of the exact potential and small distur- 
bance equations. The most widely used of these codes are those developed by Jameson 
and Caughey (refs. 1 and 2); FL022, the most reliable of these codes, computes the 
flow about an Isolated wing on a wall. FL022 has been used extensively In the design 
of advanced wxngs for transonic aircraft. Data from two such wings designed for the 
Vought A-7 fighter are compared with theoretical calculations from four different 
transonic codes to show the strengths and weaknesses of each code. The two wings were 
designed during a cooperative study Involving the Vought Corp. , the David Taylor Naval 
Ship Research and Development Center, and the Ames Research Center. During the 
cooperative study a wing design code was developed by coupling FL022 with a numerical 
optimization algorithm (ref. 3). This report will demonstrate that wing design for 
transonic aircraft by use of Isolated wing codes may lead to unsatisfactory results If 
the aircraft has a low fineness ratio fuselage. 


FL022 COMPARISON WITH EXPERIMENT 


FL022 Is a full potential, non-conservative, transonic analysis code which uses 
an explicit, square root transformation to map an Isolated wing to a "bump" on a 
plane. 

Theoretically predicted pressure distributions from FL022 disagreed significantly 
from test data at the wing design condition, 0.85 Mach number and 0.30C{^ (fig. 1). The 


wing performed satlefactorllyi but when compared with the original A-7» underpradlc- 
tlon of shock strength may have resulted in premature drag divergence. Possible 
causes of this disagreement were identified as code approximations. Isolated wing 
limitations, boundary layer analysis approximations, and test inaccuracies. 

A strip boundary layer displacement thickness computed for the design pressure 
distribution was added to the wing geometry during both design and analysis which 
Introduced an error into the theoretical calculations because the experimental span- 
wise pressures differ from the design pi.esaures which resulted in an Improper dis- 
placement thickness. In order to assess this error. FL022NM predictions were com- 
pared with the experiment (fig. 2). FL022NM combines Jameson's FL022 code with a 
Nash-McDonald boundary layer analysis. Displacement thicknesses were computed for 
each spanwlse strip and updated as the solution converged. FL022NM predictions were 
somewhat closer to the experiment than FL022 pressures, but the strong experimental 
shock was still not predicted. A comparison of FL022 and FL022NM pressures is shown 
in figure 3. 


FL028 COMPARISON WITH EXPERIMENT 


FL028 is a full potential, fully conservative, transonic analysis code developed 
by Jameson (ref. 3). The code uses the finite voltime approach and solves for the 
flow about a wing in the presence of a finite body. The body is mapped to a vertical 
slit by use of a generalized Joukowsky transfonngtlon. 

The A-7 body with a 12% thick supercritical wing was analyzed using FL028 on the 
NASA Ames CDC 7600 computer. Initial solutions were remarkably bad. Wing pressures 
were characterized by chordwise oscillations. Comparable irregularities and discon- 
tinuities were seen in the wing mapped surface coordinates at the root and tip. 

Since the correctly input wing-body geometry caused unsmooth transformed coordinates, 
it was obvious that some compromise in the geometric representation would be required. 
The body was modified because of its incompatibility with the code approximations. 

The Joukowsky transformation had difficulty mapping the A-7 body which has a nearly 
rectangular fuselage. Progressive rounding of the body comers had a smoothing 
effect on the root and tip mapped coordinates and on the resultant pressures. A 
slight lowering of the A-7 high wing was required by this body mouiflcatlon. These 
geometry changes smoothed wing pressures, but oscillations were still so bad that it 
was impossible to compare the theory with experiment. 

It seems illogical that smoothing the mapped surface coordinates had not caused 
a more pronounced improvement in the wing pressures. The code was modified to plot 
mapped coordinates at every computational station Instead of only the root and tip. 
While the earlier changes had smoothed the root and tip. intermediate span stations 
were still very irregular — particularly in the 302-40% semi-span region. While the, 
use of a circular fuselage model seemed the easiest ^Ix. wing height changes were 
tried first. It was found that even a slight lowering of the wing relative to the 
body made dramatic improvements in smoothing the mapped coordinates. The best 
results were achieved with a mid-wing position on the fuselage with slightly rounded 
comers. Adjusting input geometry to achieve a criteria of smooth mapped coordinates 
may lack mathematical rigor, but it does produce improved results. Further, it is a 
relatively easy to use criterion and is not dependent upon agreement of theory with a 
particular set of test data. 
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FL022 and early FL028 conpariaona with axpariaant ahowad batter agraamant If the 
theory was nm at a allghtly higher Mach number than the teat data. Other Investi- 
gators had reached this conclusion. The bast agreement with 0.85 Mach number experi- 
mental pressures was achieved with theoretical data run at 0.86 Mach number. It was 
hypothesised that this relatively large Mach shift might Indicate a small experimen- 
tal calibration error. In retrospect, after making comparisons with other codes, 
this reasoning appears specious. 

Once a satisfactory input geometry was achieved, analyses were run at experimen- 
tal lift conditions. The theory was run at approximately one degree higher angle-of- 
attack in order to match the test lift. Each analysis was run SO Iterations on the 
crude and medium grids, and 250 iterations on a fine grid of 160 x 16 x 28. A typi- 
cal average residual was converged to 2 x 10~^. Bach run required approximately 
1700 sec of CDC 7600 time. A few mns were executed an additional 250 iterations on 
the fine grid, doubling the run time and reducing the average residual to 0.6 x 10"^, 
while having little effect on the pressures. The tip shock was strengthened and 
moved aft very slightly, but pressures over most of the wing were unchanged. 

FL028 wing pressure predictions are compared with the experiment at 0.34 C^, the 
wing design condition (fig. 4). Some mild upper surface pressure oscillations are 
evident from the wing root to mid-span. Agreement of shock strength and position is 
very good. At O.SIS C^, agreement is equally good except near the tip where the 
experimental shock has shifted forward (fig. 5). Oil flow photographs indicate that 
there is probably some trailing-edge separation outboard of 90Z semi-span at this 
lift. Shock induced separation has spread to mid-span at the next higher lift coef- 
ficient 0.63 (fig. 6). It is interesting to note that small oscillations in experi- 
mental upper surface pressures appear to have a correspondence with theoretical 
oscillations. 

FL028 wing-body predictions show much closer agreement with the experiment than 
do FL022 isolated wing predictions. Is this difference caused by body effects or by 
differences in the two codes? FL028 was run with an isolated wing in order to inves- 
tigate this question. Isolated wing pressures are compared with test data for 
0.34 and 0.515 C^, in figures 7 and 8. Notice that pressure oscillations are no 
longer present. Both FL028 and FL022 predict similar pressure distribution for an 
Isolated wing although FL028 predicts a slightly stronger shock. It appears that the 
improved experiment-theory correlation observed with FL028 is due to the inclusion of 
body effects rather than code differences. 


FL030 COMPARISON WITH EXPERIMENT 


FL030 is a finite volume wing-body transonic analysis code developed by Caughey 
(ref. 2). This code uses body oriented coordinates instead of the Joukowsky trans- 
formation oi FL028. Geometry input for the 12Z wing was Identical to FL028 input 
except the nose was not closed and the wing was raised back to its original high 
position on the fuselage (fig. 9). A 9Z thick supercritical wing was also 
investigated. 

No geometry problems were encountered in running FL030. Most analyses were run 
200 iterations on a crude grid and 200 iterations on a moderate grid of 
80 X 12 X 16 points. One solution was run an additional 200 iterations on a fine 
grid, 160 X 24 X 32, The fine grid solution showed very little pressure change over 
the moderate grid solution except at the leading edge (fig. 10) yet required four 
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tlrntt the conputer tint (3720 ate vs 950 ssc). A pressurs splks sppvoxlMtsly two 
chord lengths aft of the wing root trailing edge was evident on all runs. While it 
is not clear that this anomaly degraded pressures on the wing, there is an obvious 
error in the code. 

Theory and experiment are compared for the 9Z thick wing-body configuxation in 
figures 11 through 13 for 0.28 C]^, 0.416 and 0.57 Agreement is generally 
good prior to separation except near the tip. Hodel manufacturing errors may account 
for pressure differences at the tip and at the leading-edge near mid-span. The 
effect of a 0.01 Mach nund>er change on FL030 wing pressures is shown in figure 14. 

The differences are concentrated in the vicinity of the shock and are largest near 
the tip. 

Comparison of theory with experiment for the 12Z thick wing is presented in 
figure IS for the design condition. Agreement is good, with the largest errors beinc 
in the predicted shock strength at 0.4 semi-span and shock position at 0.88 semi-span 
Results at 0.52 are shown in figure 16. As discussed earlier, some traillng-edge 
separation exists for this condition at the wing tip. 

Although both FL028 and FL030 show good agreement <rith the experiment* signifi- 
cant differences exist between predicted pressures (fig. 17). The wing input geom- 
etry was identical for both codes. FL030 was easier to use for the configurations 
investigated and did not require angle-of-attack adjustments to match experimental 
lift. 


WIBCO COMPARISON WITH EXPERIMENT 


UIBCO is a sm<iil perturbation transonic wing-body analysis code written by 
Boppe (ref. 4). A two-dimensional boundary layer subroutine is linked to this cede 
which adds a displacement thickness to the wing geometry and updates it at specified 
Intervals. The 122 thick supercritical wing was analyzed using WIBCO on NASA 
Langley's CYBER 175 computer. In order to save time, the wing was run on a circular 
fuselage in a mid-wing position. 

Predicted pressures are compared with the experiment at 0.34 Cj^ and 0.52 C|^ in 
figures 18 and 19. Theory was run at the test angles-of-attack. Agreement is sur- 
prisingly good considering the geometric and code approximations. 

Isolated wing and wing-body predictions are compared in figure 20. Relatively 
little difference is seen between isolated wing and wing-body predicted pressures 
except at the root. Significant differences were predicted by FL028 across the 
entire span. 


CONCLUSIONS 


Three wing-body analysis codes have bean evaluated by comparing predicted wing 
pressure distributions with experiment. All three codes show significant Improvemenl 
over the isolated wing code, FL022. FL028 was very difficult to use for the particu* 

lar fuselage analyzed. FL030 was such easier to use than FL028 but was more expen- 
sive. The simplest method evaluated, WIBCO, gave very good results without excess! v< 


4 


run tines. Unfortunately* ail codes investigated are too expensive to be used for 
extensive design perturbations. 
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